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Nomenclature for Multi-mesh Gear Dynamics Program 
Evaluation and Enhancements 




Symbol 


Description 


Variable Friction — Variables 


L/W 

V r 

v s 



Load per Unit Face Width 
Total Velocity 
Sliding Velocity 
Absolute Viscosity of Oil 
Coefficient of Friction 
Specific Gravity of the Oil 
Kinematic Viscosity 
Flash Temperature 
Inlet Temperature 


Runout — Variables 



n 

N 




Equivalent Tooth Spacing Error 
Input Displacement Error due to runout 
Current Tooth Number 
Total Number of Pinion Teeth 
Pressure Angle 


Flexible/Ring Carrier — Equation Variables 


m 5 =I 5 


*7 

™ci =I C, 

d* pi ,d rp . 

Kc -1 /K rl 
Ys r =©s 

y p; = e p; 

Yc* = 9 C ; 


*7 







/ 

To^ c 


n 


' x p* 

X 7 



'® r ?A 

P<p- 

K 

%' 4 \ 


Sun Rotational Inertia 

Ring, Planet, and Carrier Rotational Inertias 

Tooth Pair Damping 

Ring or Carrier Segment Stiffness 

Tooth Pair Stiffness 

Sun Tooth Displacement along line of action 
Planet Tooth Displacement along line of 
action 

Carrier Tooth Displacement along line of 
action 

Ring Tooth Displacement along line of action 
Torque Input 
Involute Base Radii 

Output Torque 

Total Number of Planets 

Deflections (After transformation to first order 

differential equations) 

Tooth Spacing Error 

Profile Modifications 

Tooth-pair Contact Identity Function 
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Single Mesh — Equation Variables 


L ° . 

^5p > ^Vp 
L c „ ,L 


Relative Displacement along line of action 

Static Tooth Load 

Involute Profile Modification 

Tooth Pair Stiffness 
Tooth-Pair Meshing Loads 


Subscripted Notation for Equation Variables 

r Ring 

s Sun 

c Carrier 

p Planet 

b Base 

sp Combination of Sun and Planet 

rp Combination of Ring and Planet 

i Mesh number 

j Tooth Pair Number 


I. INTRODUCTION AND SIMMARY 


A. background 

A multiple mesh gear dynamics computer program has been continually 
developed and modified during the last three years under NASA funding; see 
Reference 1. The original single mesh program was expanded to include 
epicyclic gear systems with internal, buttress, or helical tooth forms. This 
NASA contract includes the following modifications: variable contact 

friction, planet cage and ring gear rim flexibility options, user friendly 
options, dynamic side bands, speed survey option, and combining the single 
mesh and multiple mesh options into one general program. In addition to the 
program modifications, the program was evaluated by comparing calculated 
values to published test data and to test data taken on Hamilton Standard 
turboprop reduction gearbox. 

b. Analysis Approach 

The basic assumptions of References 1 and 2 for the theoretical 
derivations apply in addition to the following assumptions for the flexible 
planet carrier and ring gear rim development and the sun-planet(s) or ring- 
planet(s) options: 

1. The flexible ring gear rim and carrier flexibility option is assumed 
to be only along the 1 ines-of-action of the respective gear meshes. 

2. The static output torque is assumed to be divided evenly among 
either the planet carrier segments or the ring gear rim segments, 
depending on the type of planetary system. 

3. Any radial component of motion is assumed to he secondary to the 
tooth-pair motion along the 1 ine-of-action . 

4. The damping of carrier and ring rim segments is assumed to be 
insignificant relative to the tooth-pair meshing damping. 

5. For the non-pl anetary options, the input torque and speed are assumed 
to be at the sun gear for the sun-planet(s) option and at the 
planet(s) for the r ing-pl anet(s ) option. 

6. For the non-pl anetary options the planet carrier is assumed rigid and 
fixed or grounded. 

C . Summary 


The first program modification incorporates a more accurate, variable 
contact coefficient of friction model. This utilizes the method of Kelley and 
results in lower flash temperatures, provided the constant coefficient of 
friction previously assumed was .05 to .08. 


The next modification adds flexibility to the planet carrier and ring 
gear rim by introducing azimuthal (along the 1 ine-of-action) degrees-of- 
freedom between each planet for each carrier and ring rim section. This 
results in more than doubling the number of first-order equations to be solved 
numerically, i.e. (6n +2) equations, where n is the number of planets. The 
effect of this modification on the dynamic results has not been resolved. The 
check cases run to check the analytical model and computer code resulted in 
numerical instabilities. Cases that had been successfully run without the 
added degrees-of-freedom were difficult to obtain converged solutions for with 
reasonable carrier flexibilities. The analytical approach and system 
equations were successfully checked by obtaining average system natural 
frequencies from the eigenvalue solution and comparing against known results. 
The problem most likely is associated with the code used for the numerical 
solution of the system equations. 

Another program modification makes the program more user friendly. 

Changes were made to allow for varying levels of input. The basic level 
requires only the basic geometric and operational data (rpm, torque) with the 
other values being defaulted. At the higher levels more and more of the 
detailed information is required. This type of input should be helpful to the 
gear designer where at the early design phase the basic level (level 1) is all 
that is required. Whereas at the final design phase all three levels will be 
required. In addition, the output was refined for better understanding of the 
results . 

The dynamic side band option modification allows the user to investigate 
the effect of run-out error on the input pinion (sun) for a single mesh of two 
external gears. A sinusoidal error distribution is applied after the tooth 
error solution is obtained, to simulate runout error. The solution is 
obtained for each tooth mesh and the maximum summary tables are output. The 
loads are saved and can be used for a Fourier analysis to obtain a frequency 
spectrum for sideband evaluations. 

Another program modification was the addition of a speed survey option 
thereby allowing the user to investigate tooth pair loads for a range of input 
speeds. The maximum load is located and the corresponding speed is used for 
the remaining stress calculations. 

The final modification was to incorDorate the option for a single mesh 
into the multiple mesh proqram for both spur and helical gears. Both 
sun-planet(s) and ring planet(s) options were added for either spur or helical 
gears. The sbur gear cases were compared to the closed form solutions from 
the single mesh program with the results givinq excellent correlation. Thus, 
the validity of the multiple mesh program's numerical solution was also 
verified. The combination of the two programs makes a more complete package 
of the analyses. 
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The first test case used for the program evaluation was a Stoeckicht 2K-H 
planetary gear system with three planets and spur gear teeth. Strain gauge 
data was available throuqh published articles (JSME) for three different 
static load cases; References 6, 7 and 9. The dynamic loads predicted by the 
theory were within the reported range of + 2 x standard deviation (i.e., 95%) 
of the experimental data. The theoretical loads were near the high side of 
the experimental load range for the hiqhest static load, near the middle of 
the range for the middle load, and near the lower experimental limit for the 
lowest load case. In each case, the theory underestimates the loads at high 
speeds but this is most likely due to the sun gear's reported effective mass 
being too large, causing the peak resonance to theoretically occur at a lower 
speed. Various combinations of tooth errors were applied to the two lower 
load cases resulting in good correlation between the experimental load range 
and theoretical load range. 

In addition to the non-synchronous system, a similar planetary with 
synchronous meshing was also tested with high static load. The theory 
correlated well with the experimental range. However, the test results 
indicated that a non-tooth meshing natural frequency existed in the operating 
range which influenced the dynamic tooth loads. The analytical model is for 
tooth meshing responses only and, therefore, could not predict this mode. 

The second system used for program evaluation was a three planet 
planetary used in the PT-6 series of aircraft engines. The strain gauge data 
available was for a very limited speed range. The theory predicted stresses 
considerably hiqher than the reported experimental stresses, more so for the 
higher load case. The lower load cases had better correlation with generally 
conservative theoretical predictions. Various effective masses were tried to 
examine the effect cn the peak responses. The result of decreasing the masses 
was to shift the speeds at which resonance occurred to hiqher speeds. The 
effect resulted in better correlation between the test and the theoretical 
data. However, a basic problem exists when comparing analytical to 
experimental stresses of gear teeth for dynamic response. Strain gauge size 
and location can influence the results significantly. Frequently the gauges 
are used as the load transducers by calibration from a known static load. The 
PT-6 data, though, was reported as stress only. 

The last system compared was a Hamilton Standard first-stage turboprop 
reduction gearbox. Limited gear tooth strain gauge data was available for 
four levels of loading. For the lowest load case, the theory yielded low 
estimates of the stresses, especially at higher speeds, while for the higher 
load cases the theory predicted conservative stress levels. The strain gauge 
data had to he used directly because a static load calibration had not been 
done. Therefore, the same problems exist for this data as with the PT-6 
data. A speed survey was done comparing loads with and without profile 
modifications. The speed survey illustrated the effect of a fixed center 
input pinion versus a floating pinion. The analytical model assumes a fixed 


center and with the non-synchronous meshing resulted in an unequal load share 
between the planets (idlers) for portions of the operating range. The actual 
test setup had a floating pinion which resulted in nearly equal load share 
Therefore, the analytical capability of a floating sun should be given 
consideration for future development. 

In general, the correlation between the test data and the program was 
good. The program predicted the general location of the natural tooth meshing 
frequencies, which in itself is a major accomplishment because of the 
nonlinear meshing action of the teeth, i.e. the non-synchronous meshing, the 
nonlinear tooth-pair stiffness, and the variable contact ratio. The program 
also indicated that non-synchronous meshing in planetary systems indeed does 
minimize the dynamic response of the system because of the varying time at 
which the respective meshes come into contact. In fact, the program showed 
that with non-synchronous meshing and contact ratios greater than 1.0 the 
dynamic tooth load/static single tooth load factor could be designed to be 
always less than 1.0 because at any instant in time the input torque would be 
divided by more than the number of planets. In fact, the average number of 
teeth in contact is the contact ratio times the number of planets. 
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II. PROGRAM MODIFICATIONS 


A. Variable Contact Friction 


Prior to this contract, the multiple mesh program required a constant 
coefficient of friction to be inputted. Modifications have been made to the 
program to internally calculate the coefficient of friction as a variable 
throughout the mesh. The reauired input now is the type of oil used. 


The coefficient of friction is a function of the Benedict-Kel ley 

parameter^) ; _ 

u.W? 


x = H* 


L/W 


( 1 ) 


where L/W is the load per unit face width, Vj is the sum velocity. Vs is 
the sliding velocity, and U 0 is the absolute viscosity of the oil. 

The empirical equation of the line representing the data points presented 
by Benedict & Kelley^) and shown in Figure 1 is: 


- . 0/30SSC0 X > . 


( 2 ) 


where^# is the coefficient of friction. 

The absolute viscosity, U 0 , used for the above expression must be input 
to the program via an oil type. The other variables are already available in 
the program. The absolute viscosity is given by: 


V. = / 2 * ( 3 ) 

where is the specific qravity of the oil and Z* is the kinematic 
viscosity, Dlotted in Figures 2 and 3 vs. oil inlet temperature. 

The input parameter has been changed from a constant coefficient of 
friction to an oil type code. A code of 1 is for MIL-L-23699 oil, a code of 2 
is for M I L-L - 7808 oil , and a code of 0 wil 1 defaul t to MIL-L-23699 oil . There 
are also codes 3 and 4 which are available for two other oil types that can be 
added later. 

The coefficient of friction influences the flash temperature, Tf, by 
the equation: 


where is the inlet temperature and M is a function of several parameters 
which can be found in Reference 1. 

A comparison case was executed using a single sun-planet mesh. The 
calculated coefficients of friction are less than the .05 assumed for the 
original case, thus the new flash temperatures are lower and should be more 
realistic. Figure 4 shows the flash temperatures before and after the 
friction coefficient modifications. 

The modified program also calculates an average coefficient of friction 
throughout the mesh for each sun-planet and ring-planet mesh. 

b. Flexible Planet Carrier and Ring Gear Rim 

The previously developed dynamic load equations for planetary gear 
systems assumed rigid planet carriers and ring gear rims. The only 
flexibilities modeled in the planetary dynamic model were the tooth-pair 
flexibilities. Presently designed aircraft type planetary gear systems 
generally have flexible planet carriers and ring gear rims because of weight 
requirements. It was felt that the gear tooth dynamics would be influenced by 
these flexibilities. Therefore, the dynamic load equations and program were 
expanded to include the planet carrier and ring gear rim flexibilities. 

The equations defining the planetary model are based on relative 
tooth-pair motions actinq along the respective meshing 1 ines-of-action . For 
the effect of the planet carrier and ring gear rim flexibilities on the gear 
tooth dynamics, the motions of the planet carrier and ring gear rim were 
considered only along the meshing 1 ines-of-action . In other words, azimuthal 
motions of the carrier and ring rim were used and any radial component of 
motion was assumed to be secondary to the tooth-pair meshing flexibility. 

Figure 5a is a schematic of the dynamic model. Azimuthal springs are 
introduced between each planet to represent the carrier and ring rim 
flexibilities. Each azimuthal carrier and ring rim segment also has a 
respective inertia. 

Dampinq of the segments is assumed to be insignificant relative to the 
tooth meshing damping. Tooth meshing damping has been found to be 
approximately 5-10% of critical and ordinary structural damping is less than 
2 % of critical . 

Using the model of Figure 5 the equilibrium equations can be written, 
for the sun gear 


fil. 


* n + 2 d^. c y, - y p . - r) +£'k f ,a- y p . -Y Ci y%k 

C»/ 77 >«.* 7 7 v /Pa, 


(5) 
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( 6 ) 


(7) 


For the planet gears 

% Yp. -<*?• (Ys - Ypy ' & ) - -*VS- (?* - Y v - ty) 

-k*P/(Ys -Ye, - Yc,) rkrp; (Yp- -Yc; -Yr/) * ^ • 

For the planet carrier segments 

/Vc/ /c/ -cJ*?'- (Y±- Yp; - ?c;) - ^ r P; ( Yp; ~Yc; ~ Ky) 

-A's />• (Ys-Yp; “Yc;) - krpp ( Ypy -Yc; -Yr;) 

Yki; (Yc; - Yc,-.,) + k C/ ., (- Yc,,, + Yc,) - • 

For the ring gear segments 

Mr,- fry -drp,.(Yp,. -fry) ~ krp,. (Yp,- -Yc, -Ky) 

* k ry„ (- Yr,.,, *Yr.) +Kr, (Yr; - Yr,„) = - ■ 

The rotation is that used 1 n References 1 and 2 and defined In Figure 5. The 
second order differential equations are solved numerically because of the 
large number of equations, (3n ♦ 1), where n Is the number of planet gears. 
The routine used requires that the differential equations be transformed to 
first order differential equations. This Is done by the following 
transformation. 

Let X 5o = £ , Ypy, • Yp,- , Yc,-- Yc; , X r,- 0 • Yr; 

Xs,-Ys , X , % , Xc,-, * Vc; , Yr it * Yr; 

Xs/* Ys , Xp;, " Yp; , Xc,-t* Ycy , X ry, * Y r y (9) 

Substituting (9) into (5) throuqh (8) qives 

"Yp.-, -Xc//) k%p t . (x^ - X P/o - X 

7/4/ 


( 8 ) 




( 10 ) 


kp /t - ^ dtp; (*S - *p; t - X C/,) - JUp. d'?/ (xp; t “ Xc// “Xr;, ) 

~ m f . K ip > (X ” x *> " ) * kTp,. ‘ K ’- i ‘ ' Xr ") =C ’ ( ' 


( 1 ?) 


*-•' - k <*••-**• - x ‘k- 4- **'(**. -*» 

- - h fir's (x,„-x. c,. 

- ,*<•! .-5S. 

*<■,< - -k- dr P; -Xc/, •*'-•')-■&,■ fi'S/i*';. -*«i 

+ k-k'-' ( ~ * Xr; -) * it kr> ' (x '“ " Xri ^y-^ kr 


(13) 



*JW« **i 

(14) 


Xp;o =•*/>„ 

(15) 


*c lo * Xc. i 

(16) 

and 

X r ,c * Xr; t 

(17) 


for a total of (6n+2) first order differential equations. The output toraue 
is assumeo to be divided evenly among either the planet carrier seqments or 
the ring gear rim segments depending on which type of planetary is beinq 
modeled . 

The tooth-pair load terms of equations (10) through (13) 


are 

- x cy.) 

(18) 

and 

(x*.-x«.-x v .) 

(19) 
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Tooth spacing error, profile modification terms, and a tooth-pair contact 
identity function must be added to the tooth-pair load terms as shown 


(x». -X/& -Xc/, - £«/>• - 4 > p ; (20) 

and f\rP/ (X/% -Xt/, (2>) 

The tooth-pair contact identity functions 4*P,' and are either 1 or 0 
depending on whether the individual tooth-pairs are in contact or not. 

Contact is determined by the value in the parentheses of (20) and (21). If 
the value is positive, the tooth-pair is in contact. If the value is 
negative, the tooth pair is not in contact. 

When a planetary system is being modeled, the ring gear is grounded and 
the v terms vanish. When a star system is being modeled the planet carrier 
is grounded and the Xt; terms vanish. 

The above development results in more than doubling the required number 
of equations. The previous solution. Reference 1, used the relative 
tooth-pair displacements as the generalized coordinates, whereas this 
development uses the absolute tooth-pair displacements as the generalized 
coordinates. The absolute displacements must be used because of the 
additional degrees-of-freedom needed for the planet carrier and (or) ring gear 
rim segments . 

C . User Friendly Options 

Some "user friendly" options have oeen incorporated into the program. 

The user can choose to run at any of three levels. The first level requires 
minimal input such as diametral pitch, pressure angle, number of teeth, face 
wi dthr , torque, speed range, equivalent masses, type of system, number of 
planets, and contact ratio if it is greater than 2. The remaining variables 
default fsee user manual for values) allowing the user to obtain approximate 
results with minimum input. The next level requires the input of level 1 plus 
damping, flash temperature data, iteration data. 2 nd planet phasing 
constants. Other optional variables that may be input at level 2 include: 
profile modi f '■at ion data, face width crowninq. Hertz stress for compliance 
calculation, Ciare stress or plane strain trigger, active face widths, 
convective compliance matrix (helical gears), geometric data, geometric 
tolerances, helix errors, tooth pair spacing errors, and initial boundary 
conditions. These level 2 optional variables will default to zero, or the 
values indicated in the user quide if not specified. Level 3 includes all of 
the above, us well as profile modification shape factors, compliance 
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constants, and 3-dimensional factors. In addition to the choice of levels, 
the documentation has been expanded to ease the user's task of setting up the 
input data for any level. The program output has also been modified for 
clarity. For further input details and default values, refer to the user 
manual . 


D. Dynamic Side Bands 

Dynamic sidebands are observed in frequency spectra where they usually 
occur in plus or minus integer multiples of the primary gear mesh frequency. 
This phenomenon is often associated with inherent manufacturing errors, such 
as gear runout . 

For the sideband option, runout errors are simulated by applying a 
sinusoidal error distribution on the input pinion. The user inputs the 
displacement deviation in center location of the sun. This is then used to 
generate a tooth spacing error array using the following equation. 


Ar 5 --4, S-V (lit ( £)) sW. 


where 


and 



is the spacing error 

is the input displacement, error 

is the current tooth number 

is the total number of teeth on the pinion 

is the pressure angle 


( 22 ) 


The error array internally generated is printed after the no-error solution 
and is followed by error passes illustrating the various tooth pair meshes 
with corresponding summary tables of maximum values output. 

In order to obtain sideband information, the load data must be used in 
conjunction with a Fourier analysis in order to look at the results in the 
frequency domain. The program outputs all of the loads and corresponding 
times, in a WRITE (7,5029) format so the data can be analyzed separately. 

An example case was executed with a displacement error of .001 inch. The 
maximum normal load varied from 3045. lbs. to 5723. lbs., while the no-error 
maximum load was 4387. lbs. The load data was used in conjunction with a 
Fourier analysis program, the results of which are illustrated in Figures 6a 
and 6b . 

E. Speed Survey Option 


The program has been modified to include the option of surveying a range 
of input speeds. The user inputs the initial and final speeds of the range 
desired, and the number of intervals into which the range is to be divided. 
The proqram steps through the speed range determining the maximum tooth load 
for each speed, where step size = (final -initial speed )/number of intervals. 


10 


The difference between consecutive loads is calculated. A change in sign from 
positive to negative indicates a peak has been passed. The speed step size is 
then divided by five, and the prior two intervals are reanalyzed to find the 
peak value. The program then continues through the range checking for other 
peaks. At the end of the range, the speed corresponding to the overall 
maximum load is retained and used for the remaining calculations. Examples 
can be seen in Figure 23 and Figure 30. 

The user should note that a peak could be missed if the number of 
intervals specified is too small; however, CPU time will increase 
proportionally with increased number of intervals. A missed peak is likely to 
occur for cases with low damping, which would lead to narrow sharp peaks. One 
indication of a missed peak would be if the loads for the smaller step size 
continually decrease through the interval. If this occurs, the user should 
re-run the program with different intervals. One should also note that the 
boundary conditions are iterated for each RPM; thus, if the number of boundary 
condition iterations is large, a smaller range should be used with user 
specified boundary conditions to avoid using excessive CPU time 

There are also two messages that may appear in the printout. If the load 
monotonically decreases in the first speed interval, a message is printed that 
the initial rpm is on a decreasing load curve. Similarly, if the final speed 
is on an increasing load curve a corresponding message is printec. These 
messages indicate that the user may want to re-run the program anc check 
either lower or higher speeds, respectively. These are in addition to loads 
for each mesh at each speed. 

F . Non planet ary Systems 

1 . Sun-Planet(s) Option 

The original multiple mesh program simulated the sun-planet(s) option by 
the user inputting a trigger; SNPLBT = 1.0, which would set the compliance of 
the ring teeth to a very large value. The input now requires the system type 
trigger, K, be set to 4 for the sun-,> anet(s) option. This option is now 
solved exactly with the system equations given below. 

The sun is the input pinion while the planets are the output gears as 
shown in the following figure: 





The following equilibrium equations can be written corresponding to the 
prior figure: 


Is Ys + <£- d*?.' Ysp. + <£. A’sff’ Ytp- * - 5 — 

/•/ ' /V 7 / ^ 


(Sun gear) 


(23) 


(Planet gears) 


I P; Ypy Ysp; “tep r Ysp t * - 

where andtf/r^are the viscous damping terms, I s and Ip^ are the 
rotational inertias, A’y/ andyfcp,- are tooth pair stiffness terms, T-j n and 
T ou tpi are input and output torques, and R^s and R^ are gear involute 
base radii . 


(24) 


Subtracting (24) from (23) yields n relative motion equations: 



For a one planet system (n=l), equation (25) reduces to Richardson's 
equation! ■) for a single external -external gear mesh as follows: 






4- 


TooT 


7 TV 



~7oUT ^ 

#bp 



/ 


/»g 


(26) 

(27) 


(28) 
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where is the meshing effective mass and L 0 is the static tooth 
load consistent with the gear torques. 

For multiple planet cases equation 25 can be rewritten as: 

% * h \ ' % • =»• 

* ^ 

\f>. * Yifiy -e. 


- z 


Here and 


i 2 - 

>sp 


are the tooth spacing error and the involute 


(29) 


(30) 


(31) 


profile modification terms, respectfully, of mesh c and tooth-pair^. 

The tooth-pairwneshing loads are — z 2 ”) 

Ur; ■ <E (■ L % ^ - J ^ * ) 

where m is the number of teeth in contact at mesh i , *1 sp is the tooth pair 
stiffness and ^ is a tooth-pair contact identity function. 


(32) 


Transform to 2n first-order equations of the following form as 

was done above Tor the flexible carrier and ring gear rim equations of motion: 

LU P . + -L Up. (33) 

1 ® 


' *V,y 


(34) 


These equations are similar to the sun-planet equations of a star system. 
Reference 1, except there are no ring-planet terms or equations. 

A test case was executed for spur gears using both the single-mesh 
program, which has a closed form solution, and the multiple-mesh program, 
which uses a numerical solution, with one planet. Table 1 contains the 
solutions from the multiple mesh program compared to the single mesh solution 
and the corresponding percent differences. 
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For a boundary condition tolerance of 0.5%, the differences between the 
closed form solution and the numerical solution maximum values were less than 
2*. This comparison also verifies the validity of the numerical solution 
iteration parameters used in the multiple mesh program. Figures 7 and 8 show 
the single mesh loads and multiple mesh loads with the corresponding data 
sets, both for one mesh cycle, respectively. 

2. Ring-Planet(s) O ption 

This ring-planet is input is similar to the sun-planet option in that the 
old trigger, 5NPLBT = 2.0, would set the compliance of the sun teeth to a very 
large value. The input now required .is that the system type trigger, K, be 
set to 5 for a ring-planet(s) mesh. 

The planets act as the input pinions, with the input torque being per 
planet, while the ring acts as the output gear’as shown below: 
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The equilibrium equations corresponding to the previous illustration can 
be written for the ring gears and planet gear as folTows: 


^ t - £ +* Ye* - £ krp. Yr f -- ^ 




Y,- * Yrf; + Krp,- Yrpy 


- 


(Ring gear) 
(35) 

(Planet gears) 
(35) 


Subtracting (35) from (36) the following relative motion equations can be 
obtained: 
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Equation (37) can be rewritten as 
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Transform equation (38) to 2n first order equations: 
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These equations are similar to the ring-planet equations for a star 
system. Reference 1, with no sun-planet terms or equations. 

A test case was run on both the single mesh and multiple mesh program for 
an internal -external mesh. The maximum values are presented in Table 2, along 
with the percent differences between the exact single mesh solutions and the 
numerical solutions of the multiple mesh program. The table is followed by 
the input data set for the multiple mesh program. Figure 9. Again there is 
excellent correlation, with the maximum difference being 3 .85% . The higher 
percentage difference, compared to the sun-planet mesh, is due to a 1% 
convergence tolerance for boundary condition iterations, rather than 0.5% used 
for the sun-planet mesh. 

3 . Profile Modifications 

In addition to the options eliminating the planet carrier and ring or sun 
gear, an additional input variable has been added for tooth profile 
modifications. The user can now specify the amounts of tooth tip and flank 
reliefs for each gear for a more comprehensive tooth profile modification 
table given in the program output. 
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III. EVALUATION OF PROGRAM 


A. Introduction 


Task II required evaluating the analytical method by comparing the 
program results with published test data and Hamilton Standard test data from 
a turboprop reduction gearbox. After extensive searching through engineering 
indexes via computer, and by hand for older publications, it was found that 
very little published test data of sufficient detail existed. An exception to 
this was a series of articles published in the Bulletin of the Japanese 
Society of Mechanical Engineers. They have published several articles which 
document their experimental investigations into dynamic tooth loading in 
Dlanetary gears. 

Another article, reference 4, was located that implied strain gauge data 
had been taken at Pratt & Whitney Aircraft of Canada on a planetary system. 

The author was contacted and asked to send test data. The data received was 
very limited in that only a small speed range was tested and there was little 
information as to the specific test procedure. Limited data was also the case 
with the Hamilton Standard test program for the turboprop gearbox. The strain 
gauge data on the gear teeth was limited due to the number of available 
monitoring channels . 


B. Experimental and Modeling Characteristics 

1 . Japanese Stoeckicht 2K-H Planetary 

Extensive testing was performed on a nonsynchronous meshing, three 
planet, planetary spur gear system. The basic characteristics are found in 
Table 3 and an illustration is given in Figure 10. Table 4 contains the 
character istics of a similar test system that has synchronous meshing. 

The nonsynchronous meshing system was designed for 20KW of transmitted 
power, a speed ratio of 6, and the working range of the higher speed shaft 
from 1900 to 7200 rpm. The synchronous meshing system has similar design 
features . 

A power circulating type planetary gear testing machine was used to 
monitor the fillet strains in three sun gear teeth. Resistance wire strain 
gauges were placed on the compression side of the roots of tooth numbers 1, 5 
and 10 of the sun gear. A tooth of the sun gear will engage 30 times with the 


-5 


17 


planet gears before returning to the same engaging position. Therefore, 90 
pieces of fillet strain data, 30 per tooth, were used to obtain a mean dynamic 
load factor and coefficient of variation. 

Typical input data sets are shown in Figure 11. The tooth geometry was 
obtained from the program geometric preprocessor as this information was not 
published. The program geometric preprocessor generates the basic tooth 
geometry using basic involute relationships given in Reference 1. A model 
damping ratio of 10% was used, corresponding to the 10% ratio used in the 
Japanese theoretical calculations. The equivalent mass for the planet carrier 
includes additional mass due to a low speed rotor in the test set-up. 

2 . First-Stage Reduction gear. Model A-27 

The Pratt & Whitney epicyclic, first-stage reduction gear for the PT6 
series of engines has 3 planets with synchronous meshing. The basic 
characteristics are shown in Table 5 and the system is illustrated in Figure 
12. The input speed is 30,000 RPM or higher for normal operation! 4 ) . 

Two strain gauges were placed in a fillet of the sun gear, each located 
1/4 inch from the edge. The number of fillets monitored was not available. 
Data was taken for several different tangential load cases and for various sun 
gear rotational speeds. 

An input data set is illustrated in Figure 13. The detailed tooth 
geometry was available for all the gears; thus, the geometric preprocessor was 
bypassed. Modal damping ratios of 20% were used for the various loads, and 
two load cases were also run with 10% damping ratios. 

3 . Hamilton Standard Turboprop Gearbox 

The turboprop gearbox is a two-stage non-planetary system with 
characteristics found in Table 6 and a system diagram found in Figure 14. The 
tooth fillet strain gauge data available was for the first stage only. The 
gauges were located on an idler gear in a total of 6 fillets, 3 consecutive 
fillets on opposite sides of the gear, with 2 gauges located in the center of 
the width and 2 gauqes eacn on both edges, as shown in Figure 15. Maximum 
stresses were observed for several loading conditions at various speeds; 
however, the maximum number of these gauges observed for a certain condition 
was four and the minimum two. This was due to a limited number of monitoring 
channels. For better load share between the two idler gears a floating input 
pinion (sun gear) was used for the design. 

An example input data set can be found in Figure 16. Damping ratios of 
5% and 10% were used. The geometry was directly input, bypassing the 
preprocessor, with the actual profile modifications being included for most of 
the cases, while one case was run with no modifications for load comparison. 
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Data & Results 


1 . Stoeckicht 2K-H PI anetary 

The test data for the Stoeckicht 2K-H planetary system was presented as a 
mean dynamic load factor and a coefficient of variation. The dynamic load 
factor is defined as the ratio of the dynamic strain to the static strain at 
corresponding positions, and the coefficient of variation is the standard 
deviation divided by the mean dynamic load. These data were then used to obtain 
the dynamic load range that included +_ 2 standard deviations, or the range 
where approximately 95% of the loads would be expected to occur based on the 
experimental data. The test data used were the maximum fillet strains in a 
zone of single tooth contact. 

Figure 17 illustrates the maximum dynamic loads for a static tangential 
tooth load of 130N (29.25 lbs.) for the nonsynchronous meshing case. The 
original effective masses are the same as the authors w) used in their 
evaluation, with the carrier mass including the low speed rotor effective 
mass, as is the damping ratio of 10%. Variations of both of these parameters 
are also shown. The Japanese theoretical calculations were also presented as 
a mean dynamic load factor and a coefficient of variation, where the mean 
dynamic load factor is apDroximately 1. The upper limit of their theoretical 
variation (+2 standard deviations) appears on Figure 17 and includes the 
theoretical effect of tooth profile errors. 

Figure 18 illustrates maximum dynamic loads for a static tangential load 
of 90N (20.25 lbs.), while Figure 19 shows results for a static tangential 
load of 50N (11.25 lbs.) Both figures are for nonsynchronous meshing and 
include computer runs for both with and without tooth spacing errors. Figure 
19 also includes an effective mass variation for both with and without tooth 
spacing errors. 

The synchronous meshing case is shown in Figure 20 in which the 
conditions are similar to chose used in obtaining Figure 17. The tooth 
numbers and diametral pitch were adjusted while the pitch radii, oressure 
angle and speed ratio were keot constant. A case with a tooth error is also 
illustrated in Figure 20. 

2 . A-27 First-Stage Reduction Gears 

The test data obtained for the PT-6 first-stage gear, model A-27, 
consisted of a plot of maximum peak sun gear dynamic stress vs. sun gear 
rotational speed for various constant transmitted loads. No information was 
available as to how many teeth were strain gauged, except a diagram indicating 
only two gauges on one tooth root and reference 4 indicating two sun gear 
teeth were strain gauged. Thus, the variation of stresses was not obtained 
experimental ly. 

Figure 21 contains the lower load cases, for transmitted loads of 445N, 

890N and 1335N (100 lbs., 200 lbs., and 300 lbs.). All three cases were run 
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with a relatively high modal damping ratio of 20% and the 1335N (300 lbs.) 
load was also run at 10% damping. Similarly, figure 22 illustrates the higher 
tangential load range. Loads of 1557N, 1668N, 1779N and 2002N (350 lbs., 

375 lbs., 400 lbs. and 450 lbs.) were all run with 20% damping, and the 2002N 
case was also run with 10% damping. 

Figure 23 shows maximum normal loads vs. speeds for several effective 
masses and for both sun-planet and ring-planet meshes. All load curves are 
for the 2002N (450 lbs.) tangential load case. The sun-planet mesh is shown 
to have two resonant speeds, while the ring-planet mesh has one major peak 
that occurs at a speed in between the two sun-planet peaks. 

3 . Hamilton Standard Turboprop First-Stage Reduction Gear 

The strain gauge data taken on the turboprop reduction gear system was 
limited in applicability for comparison to the program results. One idler 
gear, which meshes with the input pinion had a total of 6 gauges in the tooth 
root area. For each run at various torque levels, different combinations of 
gauges and corresponding peak stresses were monitored. Gauges had been 
installed on the bull gear teeth also, but they were not monitored due to the 
lack of channels available. Thus, no data was available for the second-stage, 
only for the first-stage idler gear that meshes with the input pinion. 

Figures 24 and 25 show the average of the experimental peak stresses and 
the average of the two theoretical input pinion-idler meshes. Cases were run 
with profile modifications, as the actual gears tested had profile 
modifications. The effects of different damping ratios are also illustrated. 

Figures 26 through 29 show each load level individually, with maximum, 
minimum and average experimental peak stress points plotted. Different planet 
(idler) stresses calculated by the multiple mesh program, due to a fixed input 
shaft in the model and nonsynchronous meshing, are the boundaries for one 
theoretical band. Note: as discussed earlier, the program assumes fixed 

centers for the gears, whereas, the input pinion of the system is floating 
between the idlers. The other theoretical band arises from a nonuniform load 
distribution across the face width due to the location of the web. By use of 
finite element analysis done at Hamilton, the loads were found to be 
approximately 12% greater at the center of the face. 

Figure 30 is a speed survey of idler gear load vs. input pinion speed for 
the 30% torque case which shows two speeds where resonance occurs. This also 
illustrates the unequal loads that can result due to the theoretical model of 
a fixed, rather than floating pinion. The differences between loads on a 
standard involute form and a tooth with a modified involute can also be seen. 
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0. Discussion of Results 


1 . Japanese Stoeckicht 2K-H 

The Japanese authors performed extensive tests on the Stoeckicht 
planetary gear system. Two systems were tested, one synchronous meshing and 
one ncnsynchronous meshing. The two systems were kept as similar as possible 
with the major changes being diametral pitch and number of teeth, in order to 
make a direct comparison oetween systems with the major difference being the 
planet phasing^) . 

Figure 17 shows the nonsynchronous meshing system which parallels Figure 
20 illustrating the synchronous meshing system. The range of experimental 
strain gauge results indicate a resonance near 2600 RPM, which is due to a 
resonance of the total torsional system resulting in torque variations (6 ) . 

The analytical program assumes constant torque, thus this resonance would not 
he predicted. The next resonance indicated experimentally is in the 5000 RPM 
vicinity. The theory predicts a resonant speed at 5200 RPM. 

There is also an experimental resonance at approximately 7200 RPM, which 
for the synchronous case, is a multiple of the meshing frequency^). The 
experimental range in Figure 20 does not indicate this significantly; however, 
the in-phase theory predicts peak loads near this speed. Synchronous meshing 
will generally cause higher loads due to simultaneous meshing of the planets. 
This is shown by comparing Figures 17 and 20 where both the experimental and 
analytical loads are high for the synchronous meshing case at the resonant 
speeds . 

Variation of the analytical effective mass has the effect of shifting the 
speed at wh’ch peak loads occur. The masses used for the majority of the 
theoretical calculations are the masses calculated for the Japanese 
theoretical eval uation(5 ) . The theoretical model of Reference 5 is very 
similar to that used here. The primary differences are that the Reference 5 
model includes additional error terms and mass terms for rotor inertias but 
apparently did not include nonlinear tooth-pair meshing characteristics. 

Note: the selection of effective masses for the analytical model is not a 

trivial task in that it is difficult to decide on the amount of attached mass 
that should be included. These attached masses could be rotors, shafting, 
mounts, etc. A complete system dynamic analysis would he required to define 
the equivalent masses, which are frequency dependent. The planetary gear 
tested contained a high and low speed rotor which were modelled by adding the 
additional mass to the ring and planet carriers, respectively. However, the 
ring is fixed for a planetary system and, therefore, has zero rotational 
moment of inetia and correspondingly, zero effective mass. The actual system 
has a radially floating ring as well as a floating sun, neither of which is 
included in the theoretical model. 
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The shift in the speed at which the peak load occurs, due to altered 
effective mass, is illustrated in Figures 17 and 19. Figure 17 shows the 
theoretical peak shifting from about 5300 RPM to approximately 7400 RPM when 
all masses are divided in half. The half mass results seem to correlate with 
the experimental data better; however, there is insufficient information to 
verify either quantity as the actual effective mass. Figure 19 illustrates a 
lower transmitted load and correspondingly, the peaks are less pronounced. 

The flatter peaks make the distinction between the peaks less obvious, 
expecially in the cases with tooth errors where both bands have excellent 
correlation with the experimental range. The half effective masses seem to 
correlate with the experimental data better at speeds above 6000 RPM, 
primarily due to a shift in peak location similar to the one observed in 
Figure 17. 

In addition to the case where all effective masses are halved, a case 
where only the planet carrier mass was decreased is also shown in Figure 17. 
There is a slight increase in loads and a slight shift in the main peak toward 
a higher speed. The differences in loads are less than 3%. Thus, in this 
case, the carrier effective mass is not the primary mass of the system even 
though it is the largest effective mass of the system. 

Figure 17 also contains a few points at a much lower damping ratio of 
5*. This increases the load at 5000 RPM by 26. 5% and emphasizes the dynamic 
peak obtained near 5000 RPM. The Japanese theoretical calculations used a 10% 
damping ratio which apparently correlated well to the experimental test data. 

The effects of tooth errors on the theoretical loads are illustrated in 
Figures 18 to 20. The loads plotted for the error cases are the minimum and 
maximum, although the maximum normally occurs in only one mesh and the average 
increase in loads for the error cases would be considerably lower. The 
manufacturing errors were estimated in reference 7. The pitch error, profile 
error and lead error were incorporated in a statisti cal tota l for a sun-planet 
mesh of magnitude .003388 mm (.00013338 in.) using from 

equation (14), Reference 8, where the 's are the various types of tooth 
errors. ^ 

Figure 18 shows the theoretical results of having one tooth error in each 
or the three sun-planet meshes. At low speeds, the maximum loads predicted 
are less than 5% higher than the no error case. In the vicinity of the 
predicted resonance, the deviation increases to about 13% higher and 6% lower 
than the no error case as one would expect. The deviation decreases at speeds 
higher than resonance. 

Figure 20 shows the synchronous meshing case and has the same magnitude 
of error applied to one tooth in all three sun-planet meshes, similar to the 
cases of Figure 18. The higher tangential load tends to decrease the effects 
of tooth spacing errors on dynamic loads, especially at higher speeds. 
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Figure 19 has the result of a different configuration of tooth errors 
plotted. The errors are of the same magnitude as in Figures 18 and 20, but 
they occur on five consecutive teeth and only one sun-planet mesh. In this 
case, which is for a very low transmitted load, the increase in load is fairly 
constant over the speed range, especially for the original mass. The half 
mass case shows more variation than the initial mass case throughout the speed 
range and seems to follow the experimental range better, except between 5000 
and 7000 RPM where the theoretical band is considerably narrower than the 
experimental . 

When the experimental and theoretical load ranges for the no tooth error 
and initial mass cases in Figures 17 through 19 are compared, it can be seen 
that for the highest static load case of 130N, the theoretical loads are near 
the upper boundary of the experimental range, for the 90N case they are more 
toward the center of the range, and for the lowest load of 50N, the 
theoretical loads are near the lower boundary. 

2 . PT6 First-stage, Model A-27 

The strain gauge data for the first-stage A-27 was obtained from Pratt & 
Whitney of Canada in the form of a plot of stress vs. speed for several 
transmitted loads. This data is reproduced on Figures 21 and 22. 

The results plotted in Figure 21 are for 20% damning and predict stresses 
higher than tne experimental data for the 890N (200 lbs.) and 1335N (300 
lbs.). The 890N theoretical case predicts stresses from 17% to 58% higher at 
13,800 RPM and 21,800 RPM, respectively. The 1335N theoretical case predicts 
stresses from 28% to 139% higher at 28,600 and 22,500 RPM. The lowest load, 
445N (100 lb.) has the best magnitude correlation, especially at lower speeds 
-- up to about 23,000 RPM the difference is less than 10%. However, at 30,000 
RPM the theoretical stress is 58% lower the experimental stress. The 10% 
damping case emphasizes any peaks and thus increases the maximum differences 
and decreases the minimum differences to 155% and 21% for the 1335N (300 lbs.) 
case. Strain gauqe size and placement can significantly affect the results. 

No static tooth load calibration was made on the qauges; therefore, the qauge 
data can only be used to indicate tooth load trends. In fact, if the 
experimental data is normalized, then the correlation is qood with theory. 

Figure 22 illustrates the higher transmitted load cases. The 
experimental data indicates a resonance occurinq between 28,000 and 
31,000 RPM. The minimum deviations for the 20% damping cases range from 27% 
to 84% corresponding to the loads ranqing from 1557N (350 lbs.) to 2002N 
(450 lbs.). Similarly, the maximum deviations range from 184% to 66% for 
loads of 1557N to 2002N. The 10% damping case indicates smaller deviations of 
about 51% near the experimental resonance, however the theoretical peak occurs 
at a lower speed at which insufficient test data is available. 
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The experimental data is strange in that the maximum stresses vary by 
less than 9% near resonance, while the applied tangential load increases by 
29%. It may be that the damping of the system is load dependent. The theory 
predicts stresses that increase in a similar manner to the increase in applied 
load, with the overall stress increasing by approximately 29%, for assumed 
constant damping. 

The resonance indicated by the experimental data near 29,000 RPM could 
have been due to the influence of the ring-planet mesh or a system resonance. 
The theory predicts resonances for tooth-pair meshing, but is not set up to 
predict overall system responses or shaft responses. The sharp load peak for 
the ring-planet mesh shown in Figure 23 indicates a strong ring-planet 
resonance that would be likely to be picked up in experimental test results. 

As there is no strain gauge data for the ring or elsewhere, it is impossible 
to differentiate where in the system the resonance is occuring. 

The theoretical stresses on Figures 21 and 22 are for standard spur gear 
teeth with no profile modifications and include no tooth spacing errors. 
Experimentally, one possible reason for the measured stresses being nearly 
independent of the applied loads could be due to the effects of inherent 
manufacturing gear tooth errors. 

Figure 23 vividly illustrates the importance of the calculated effective 
mass used for the analytical model and the effect of varying the effective 
masses. As the masses are decreased, the peak responses shift to higher 
speeds. There are two resonant speeds for the sun-planet mesh predicted 
theoretically; however, experimentally only a small range of speeds was tested 
which includes only one resonance between 28,000 and 31,000 RPM. This could 
correspond to either of the peaks theoretically predicted for the sun-planet 
mesh or could correspond to the one resonance predicted for the ring-planet 
mesh . 


Figure 23 also shows the result of a lower modal damping ratio. As the 
damping ratio is reduced, the peaks become much more pronounced. One side 
effect of reduced damping ratios in running the program is generally a larger 
number of convergence iterations for boundary conditions. Therefore, in most 
cases, lower damping ratios will require more CPU time than higher damping 
ratios . 

3 . Hamilton Standard Turboprop Reduction Gearbox 

The test program for the turboprop gearbox involved strain gauginq the 
gear rims, webs, and tooth roots. At different torque levels, various 
combinations of gauges were monitored. Thus, the peak stress data from some 
of the gauges located in the fillets, see Figure 15 for locations, were 
compared to the maximum stress results from the multiple-mesh program. 
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Figures 24 and 25 illustrate the average of the experimental data points 
for each speed and the averages of the two theoretical idler (planet) 
stresses. The experimental data was averaged to eliminate the data scatter. 

The experimental data scatter is a result of gear runout, tooth spacing 
errors, tooth geometric tolerances, and strain gauge placement. Some of this 
data was actively available that could have been used but contract funding was 
limited and, therefore, beyond the scope of the study. The theoretical data 
was averaged between the two idlers (planets) because of unequal load share. 
The actual experimental setup had a vertically floating input pinion which 
resulted in equal idler load share in the test setup. 

Based on the averages, the two lower load cases. Figure 24, have good 
correlation at lower speeds, and reasonable correlation at the higher speeds. 
For the two higher load cases in Figure 25, theoretically predicted stresses 
are considerably higher than the experimental stresses for the lower speeds. 
The correlation at higher speeds is excellent with an 8.3% difference between 
the 90% torque case using 5% damping at 21,100 RPM. The 68% torque case had 
only two speeds at which experimental data was taken, thus comparison may be 
based on the trends of the other load cases. 

The theoretical data shows an increase in stress for speeds higher than 
the tested range with a peak appearing at about 24,000 RPM. If the 
theoretical effective masses were increased, shifting the peak to a lower 
speed, more of the experimental ranges would have paralleled the theoretical 
stress bands. The higher torque cases. Figures 28 and 29, show the theory 
being conservative and predicting stresses higher than the experimental data, 
especially for the lower speeds. If the peaks were shifted, the differences 
would again be decreased, yielding better correlation between test and theory. 

The effect of modal damping is shown in Figure 24. For the 50% torque 
case, the stresses vary by about 10% for an increase in the damping ratio from 
5% to 10%. The largest difference occurs as the pinion speed approaches 
resonance, as would be expected. 

Figures 26 through 29 show the minimum and maximum experimental data 
points, as well as the average experimental stress and theoretical ranges. 
However, there are extremely large variations in the experimental stresses, 
especially for the 30% and 50% torque cases which also had more experimental 
data points. As mentioned, much of the variation can be attributed to 
gearmount , tooth spacing errors, tooth tolerances and strain gauge placement. 
At the lower loads, all of these effects excluding gauge placement are much 
more influential that at the higher loads because the qear tooth deflections 
are smaller for lower loads. The tooth errors and tolerance differences are a 
larger portion of the tooth deflections. 

One theoretical band on Figures 26 to 29 represents the stress 
distribution across the tooth face due to a web located at the center of the 
face. Studies were conducted using finite element analyses to determine the 
effect the web had on the load distribution. The result of that study 
indicated a 12% increase in the center of the face and a 12% decrease near the 
edges. This variation was applied to the theoretical average to obtain the 
dashed band on Figures 2C to 29. 



Figure 30 is a theoretical speed survey for the lowest torque case, 
illustrating the effect of profile modifications. The modifications have 
little effect on the maximum loads at a resonant speed. However, the minimum 
loads prior to the resonant peak are larger for the modified profile than for 
the standard involute form. 

Theoretical unequal load sharing is also illustrated in Figure 30. In 
the operating range, the loads are fairly equal. At resonant speeds, however, 
the loads vary by as much as + 26% . The variation is less than + 6% for 
speeds of 25,000 RPM and above. The variation is due to the fixec! pinion 
mathematical model, which does not allow for a floating pinion present in the 
actual system. 
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IV. CONCLUSIONS 


1. The variable coefficient of friction modifications lead to more 
accurately predicted flash temperatures. 

2. Modifications to make the program user friendly assist the user in 
setting up the input data via three levels of input, and they aid in 
interpreting the output. 

3. The dynamic side band option allows the user to investigate the 
simulated effect of a runout error on the input pinion of a single 
external -external mesh. The number of error solutions calculated is 
equivalent to the number of teeth on the input pinion, simulating 
runout for a full mesh cycle. The load data can be used to obtain 
the frequency spectrum to evaluate the sideband phenomenon. 

4. A speed survey option allows the user to obtain load variations f^om 
speed range and to determine more readily the speeds where resonances 
are likely to occur. 

5. The addition of nonplanetary, sun-planet(s) and ring-pl anet(s) , 
systems not only allows for single-mesh equivalent systems, but also 
allows for multiple planet systems with no planet carriers, and 
results in one universal analysis program. 

6. Comparison of spur qear single mesh configurations in the multiple 
mesh program to the single mesh program, which has a closed form 
solution, indicates the numerical solution used for solving the 
differential equations in the multiple mesh program is accurate. 

7. The flexible planet carrier or ring gear rim options resulted in some 
numerical instabilities. The system equations have been verified by 
the eigenvalue solution so the problem most likely exists in the 
numerical solution routines. 

8. The comparison of the Stoeckicht 2K-H planetary system data with 
theory showed the loads correlated well. The other two systems 
evaluated had stress data rather than load data, which did not 
correlate as well with theory. In addition, the Stoeckicht load data 
was presented statistically rather than with absolute measurements 
and, therefore, corresponded more realistically to theoretical 

resul ts . 

9. The program predicted the general location of the natural frequencies 
of the three systems, even though the meshing action of the teeth is 
nonl inear . 
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10. All three test system showed a better correlation to theory at 
higher applied static load levels, while the program tended to 
underestimate the dynamic loads or stresses for lower static load 
levels. This is a result of the relative displacements of the tooth 
pair under load compared to tooth spacing errors. At low loads the 
spacing errors can significantly affect the dynamic loads. 

Therefore, for lightly loaded gears the tooth-spacing errors should 
be carefully defined and analyzed. 

11. All three test systems were shown to be very dependent on the 

equivalent masses ana damping ratio input, with respect to resonant 
speeds and load magnitudes, respectively. This phenomenon is typical 
of lowly damped systems to vibratory loads. Usually a range of 
effective mass is studied for sensitivity of dynamic response. 
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Backlash Tolerance (mm) 


TABLE 4: STOECKJCHT 2K-H SYNCHRONOUS PLANETARY CHARACTERISTICS 
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Tooth Profile error 
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HAMILTON STANDARD FIRST STAGE TURBOPROP GEARBOX CHARACTERISTICS 
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Tooth Thickness, nominal (nm) 3.90 3.35 



SPECIFIC GRAVITY 


ORIGINAL FAQ- ;; 
OF POOR QUALITY 


Figure i. Specific Gravity vs. Temperature 



TEMPERATURE , °F 



KINEMATIC VISCOSITY - CENTISTOKES 


original page is 

OF. POOR QUALITY 


Figure 3. Kinematic Viscosity vs. Temperature 



TEMPERATURE , *F 
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FLASH TEMPERRTUR FLASH TEMPERATUR 


Figure 4. Flash Temperature, Before and After Variable Friction 
Modifications to Single Mesh Program 


0 .6000E+03T 
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Figure 5. Schematic of Dynamic bedel for 
Flexible Carrier and Ring Rim 
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VALUE OF FUNCTION 


Figure 6a: Runout Loads vs. Angular Location on Pinion 

14 Tooth Planet 


THIS IS THE FINAL LOAD CASE. S4 HARMONICS 



MA6NITIJDE OF TRANSFORMED LOADS 


ORIGINAL PAGE IS' 
OF POOR QUALITY 


Figure 6b: Fourier Transformed Loads vs. Order of Harmonics 



MULTIPLES OF MESH FREQUENCY 


Mesh Frequency ^ 14 
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Single Mesh Loads and Input Data Set 









TEST 

CASE FC» 

SL’N AN!) PLANET 




3 

1 . 

2. 

8.4667 

22.5 


5 

6. 

14. 

23. 

0 . 

1.1C11 

3 

12. 

1 . 

4. 

300. 


3 

15. 

3000. 

3000. 

1 . 


4 

28. 

. 05el67 

0 . 

0 . 

.15507 

2 

52. 

. .1 

.1 



3 

140. 

25. 

180. 

.0528 


2 

150. 

005 

20. 




601. 

1 . 





0 - 1 . 
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Figure 9 


Input Data for Multiple-Mesh, Ring-Planet Option 


.TEST. 

5 

FOR RING- 
1 . 

•PLANET .OPTION 
2. 

6.4667 

22.5 

0. 0. 

5 

6. 

0. 

28. 

70. 

0. 1.1811 

4 

11. 

1.1611 

1. 

5. 

70. 

3 

_15. 

3000. 

3000, 

„_1. 


5 

30. 

12. 

.15507 



2 

52. 

.1 

.1 



1 

120. 

0. 




_ 3_ . 

2 

.... _14Q,. . 
150. 

25. 

.01 

180.. _ 

20. 

_ .0528. 

— — — _ 

1 

167. 

.0075 




1 

501. 

.0000105 




1 

541 

-.00330 




1 

699. 

0. 




0-1 
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Figure 10. Stoeckicht 2K-H Planetary System Diagram 



1 sun gear 4 coupling retainer 

2 planet gear 5 high-speed gear coupling 

3 ring gear 6 planet carrier 


Reference 7 



Figure 11. Input Data Set for Stoeckicht 2K-H System 


STOECKICHT GEAR. 






5 

1. 

2. 

8.4667 

22.5 

0. 

0. 

5 

6. 

14. 

28. 

70. 

1.1811 

1.1811 

4 

11. 

1.4184 

1 

• 

1. 

25.8 


3 

15. 

1000, 

6000. 

5. 



5 

28. 

.056167 

28.8987 

0. 

.15507 

.15507 

1 

33. 

.15507 





2 

52. 

.1 

.1 




1 _ 

.120 v _ 

0 . ... 





3 

140. 

25. 

180. 

.0528 



2 

150. 

.01 

20. 




1 

167. 

.0075 





3 

421. 

0. 

.,6666667 

.3333333 



3 

441. 

.5 

.1666667 

.8333333 



1 

699. 

0. 






0-1 
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Figure 12. A-27 System Diagram 



"V- f .. r.c 2 4 
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Figure 13 


A-27 Input Data Set 


FIRST STAGE ,_PWtC 


1 . 

3. 

19.2 

6. 

27. 

45. 

IX. 

.945 

3. 

1 *. 

949.2525 

36000.. 

28. 

.001567 

.004317 

33. 

.001017 


52. 

.20 

.20 

.120 »_ 

1. 

_ .6473 

125. 

1.18315 

2.993 

130. 

.0223 

.022 

140. 

25. 

180. 

150. _ 

.02 

30«_ 


22.5 


117. 

1 . 

1.222 

1.222 

_36000 

. 1 


0 . 

.001017 

.001017 


_ L-097 

-3.117- 

7592 

.0651 

.074 

.083 

.016 



.0528 
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Figure 16: Gearbox Input Data Set 



GEARBOX, 

INPIT PINION 

ft GEAR K 




3 

1 . 

2. 

11.0099 

25. 



5 

6. 

24. 

77. 

0 . 

2.220 

1.90 

2 

12. 

2 . 

4. 




4 

14. 

1682. 

23000. 

23000. 

1. 


5 

28. 

.0033804 

0 . 

0 . 

.017106 

.017106 

1 

52. 

.10 





4 

60. 

75.4 

71.5 

.0006 

.0006 


3 _ 

. ... 66 . . 

—15.35 

_ 84.65 

,0002 



1 

78. 

.0007 





5 

120. 

1 . 

.98515 

3.36440 

0 . 

1.19225 

4 

125. 

3.57620 

0 . 

.1534 

.1319 


2 

130. 

s_041 

.037 

0 . 



3 

140. 

20. 

120. 

.0528 



2 

150. 

.005 

20. 




2 

421. 

0 . 

.18838665 





0 - 1 . 




FIGURE 18. STOECKICHT 2K-H PLANETS NONSYNCHRONOUS 
STATIC LOAD « 90 N (20.25 LBS ) 





S9T 'oven -IVWyON 


TANGENTIAL LOAD: 



FIGURE 21 : MODEL A-27, STRESS VS SPEED LOWER TANGENTIAL LOADS 




TANGENTIAL LOAD: 



2002N(450 LB ) TANGENTIAL LOAD 



FIGURE 23: MODEL A-27, SPEED SURVEY 
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FIGURE 24. HAMILTON STANDARD GEARBOX FIRST STAGE IDLER STRESS 


RUN 43 — 90>i TORQUE 

• experimental AVG 



62 


FIGURE 25. HAMILTON STANDARD GEARBOX FIRST STAGE 



FIGURE 26. HAMILTON STANDARD GEARBOX. IDLER STRESSES-30 % TORQUE CASE 


. STRESS 
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+ E XPERIMENT AL A EXPERIMENTAL STRESS AWG 

MINIMUM AND CALCULATED AVG *12% 



FIGURE 28. HAMILTON STANDARD IDLER STRESSES— 68% TORQUE CASE 


'LA*. El 

ess. m 

FILE. ! 



INPUT TORQUE - 1 90. N • m(l 682 IN. LB) ~ 30% □ INCLUDING PROFILE MODIFICATIONS 

10 % DAMPING 

NO MODIFICATIONS. INDIVIDUAL 

NO PROFILE MODIFICATIONS. AVG. LOAD PLANET LOADS 



FIGURE 30. HAMILTON STANDARD GEARBOX, SPEED SURVEY OF LOADS 
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APPENDIX A: USER'S MANUAL OF POOR QUALITY 


*********************************** 
*********** F17S ****************** 
*********************************** 

* * USER’S MANUAL FOR MULTIPLE MESH * 

* * SPUR AND HELICAL GEAR SYSTEMS*** 

*********************************** 
********** AUG 1984 *************** 

*********************************** 

********* PROGRAM RESTRICTIONS ************** 

1. MAXIMUM PLANETS FOR SPUR AND SINGLE HELICAL GEARING 
IS 20 

2. MAXIMUM PLANETS FOR DOUBLE HELICAL GEARING IS 10 

3. MAXIMUM ITERATIONS FOP. SOLUTION CONVERGENCE IS 20 

4. MAXIMUM TEETH WITH TOOTH SPACING ERRORS IS 5 

5. 10 TEETH ARE CHECKED FOR DYNAMICS DURING TOOTH ERROR PASS 
G. MAXIMUM INVOLUTE CONTACT RATIO IS 2.0 

7. MAXIMUM NUMBER OF SUN TEETH IS GO TOR RUNOUT OPTION 

*********************************** 

* * * INPUT DATA FOR MULTIPLE MESH SPUR AND HELICAL GEAR PROGRAM * * * 

NOTE: ALL GEOMETRIC INPUT DATA IS IN TIIE ROTATIONAL PLANE AND 
ALL DEFAULT VALUES ARE 0.0 UNLESS OTHERWISE NOTED. 

ANGLES ARE IN DECREES 
FORCES ARE IN POUNDS (/•') 

LENGTHS ARE IN INCHES (IN) 

MASSES ARE IN (#-SEC**2/IN ) 

STRESSES ARE IN (///IN* *2) 

TEMPERATURES ARE IN DEGREES FAHRENHEIT 


THE DATA SET REQUIRED TO Rl’N THE MULTIPLE MESH PROGRAM IS 
DESCRIBED BELOV. AND AN EXAMPLE DATA SET WITH CORRESPONDING 
EXAMPLE PROBLEM ARE IN APPENDIX C. A LIST OF GEAR SYSTEM 
PARAMETERS WITH DESCRIPTIONS, LOCATION NUMBERS, AND SOME 
NONZERO DEFAULT VALUES FOLLOWS THE INPUT DESCRIPTION. 

CARD 1: TITLE CARD— CONTAINS CASE TITLE AND/OR DESCRIPTIVE 

INFORMATION. AT LEAST ONE COLUMN FROM 
1 TO 40 MUST BE NONBLANK. 

DATA CARDS: CARDS CONTAIN DATA IN THE FORM SHOWN IN FIGURE 1, 

NUMBER OF DATA ITEMS ON CARD, LOCATION NUMBER OF 
FIRST DATA ITEM IN LINE (SEE GEAR SYSTEM PARAMETERS 
SECTION FOR PARAMETER DESCRIPTION AM) C ORRESPOND I N(. 
LOCATION NUMBERS), AND DATA ITEMS ARE IN G FIELDS 
OF 12 SPACES. 
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PRECEDING TAGE BLANK NOT FILMED 


CASL TERMINATION CARD: TO TERMINATE A CASE THE LAST LINE MIST 

CONTAIN 0-1. IN COLUMNS 1-4. 

SUBSEQUENT CASES: ANOTHER DATA SET OF THE SAME FORMAT— TITLE 

CARD, DATA CARD, AND CASE TERMINATION CARD- 
MAY FOLLOW THE CASE TERMINATION CARD. 

JOB TERMINATION CARD: AFTER THE LAST SUBSEQUENT CASE TERMINA- 

TION CARD, A BLANK CARD (LINE) MUST BE 
INCLUDED. 

** NOTE ** ALL NUMBERS MUST BE REAL EXCEPT THE NUMBER 0* 
.ITEMS, THAT IS, A DECIMAL POINT MUST BE INCLUDED. IT IS NOT 
NECESSARY TO INPUT ZERO VALUES AS THEY WILL DEFAULT TO ZERO 
UNLESS OTHERWISE SPECIFIED. 


A BRIEF DESCRIPTION OF OUTPUT FOLLOWS THE GEAR SYSTEM 
PARAMETER SECTION. 

FOR EXAMPLE, IN APPENDIX C THE TITLE CARD READS ’ EXAMPLE DATA 
SET’. THE NEXT LINE IN THE DATA SET CONTAINS 5 PIECES OF DMA, 
STARTING AT LOCATION Itl . , FOLLOWED BY THE LEVEL (2.), THE DIAME- 
TRAL PITCH (3.4CC7), THE PRESSURE ANGLE (22.5), COAST SIDE PRES- 
SURE ANCLE AND HELIX ANCLE (BOTH - 0. FOR STANDARD SPl'!'. GEARS, 
THESE WOULD DEFAULT TO 0 IF LEFT BLANK AND THE NUMBER OF ITEMS IN 
THE LINE WERE CHANGED TO 3.) 

M0S1 OF THE INPUT VARIABLES HAVE SUFFICIENT EXPLANATION IN 
THE GEAR SYSTEM PARAMETERS SECTION, HOWEVER THE FOLLOWING 
PARAMETERS SHOULD BE CALCULATED AS FOLLOWS. 

*** EQUIVALENT MASSES: 

SUN EQUIVALENT MASS - J / (RASE RADIUS OF SUN) **2 

WHERE J IS THE MASS MOMENT OF INERTIA FOR THE SUN. THE OTHER 
COMPONENTS ARE CALCULATED SIMILARLY, USING THE CORRESPONDING 
MOMENTS AND BASE RADII. 

*** PHASING CONSTANTS: 

FOR EQUALLY SPACED PLANETS THE SUN PHASING CONSTANTS ARE 
DETERMINED BY ASSUMING THE FIRST PLANET MESH HAS A PHASING 
CONSTANT OF ZERO. THE REMAINING SUN-PLANET PHASING CONSTANTS 
ARE DETERMINED BY: 

(PLANET // - 1 )*THE FRACTIONAL REMAINDER FROM DIVIDINC THE 

NUMBER OF SUN TEETH BY Till: TOTAL NUMBER OF PLANETS 
WHERE THE 'FRACTIONAL REMAINDER’ INDICATES THE SPACING DIFFER- 
ENCE BETWEEN PLANETS. 
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THE PHASING CONSTANTS FO!! THE RINK-PLANET MESHES ARE DETERMINED 


THE SAME WAV IF THE PLANET HAS AN GDI) NUMBER Or TEETH, US INC THE 
NUMBER OF KING TEETH. IF THE PLANETS HAVE AN EVEN NUMBER OF 
TEETH THE CONSTANTS ARE CALCULATED FOR THE SUN MESHES THEN 0.5 
IS ADDED TO EACH TO OBTAIN THE RING-PLANET MESH PHASING 
CONSTANTS <KRP), DUE TO THE RING GEAR BEING .5 OFFSET FROM THE 
SUN GEAR. 


EXAMPLE CORRESPONDING TO APPENDIX C: 

SUN PLANET MESH— 14/3 - 12 2/3 

PLANET 1, KSP(.l) - 0. 

PLANET 2, (2-1) *(2/3) - .6666GG7 
KSP( 2 ) - . 6GGG6G7 

PLANET 3, ( 3- 1 ) * < 2/3 )— 1 .3333333 
KSP(3) - .3333333 


RING PLANET MESH— (EVEN NUMBER OF PLANET TEETH) 


PLANET 

1, 

KRP(l) - 

KSP(l) ♦ .5 - .5 

PLANET 

2. 

KSP( 2) 

♦ .5 - 1 . 1GGGGG7 



KRP<2) - 

. 1GGGGG7 

PLANET 

3, 

KSP( 3 ) 

* .5 - .8333333 



KRP( 3 ) - 

.8333333 


*** BOUNDARY CONDITIONS: 

AN INITIAL ESTIMATE FOR DISPLACEMENT BOUNDARY CONDITIONS MAY 
BE OBTAINED BY DIVIDING THE STATIC TOOTH LOAD AT THE PITCH RADIUS 
BY THE TOOTH SPRING RATE AT THE PITCH RADIUS. 


***** LEVELS OF INPUT 
LEVEL I INPUT 

THE FIRST INPUT ITEM WILL BE THE LEVEL DESIRED, FOLLOWED BY ITEMS 
2 THROUGH 51 AND WHERE NOTED. THE OTHER VALUES WILL DEFAULT. 

THE MATERIAL PROPERTIES, LOCATIONS 22 TO 27, DEFAULT FOR STEE.I. . 
THE TOLERANCES ARE SET TO ZERO OR DEFAULT VALUES BELOV. , ERRORS 
ARE SET TO 0. THERE ARE NO PROFILE MODIFICATIONS. OTHER NONZERO 
DEFAULT VALUES ARE: SURFACE ROUGHNESS - 25, OIL TEMPERATURE - 

ISO F, MATERIAL CONSTANT - 0.0523, DAMPING - .02, AND OIL TYPE 
MII.-L-23GDU . ANY DESIRED PLOTS CAN BE OBTAINED, SEE LOG GU1-G0S, 
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AM) A CHECK ON INPUT DATA CAN BE MADE, It 69?. IF CONTACT RATIO 
IS HIGH (GREATER THAN 2) IT MIST BE INPUT IN LOC 54 & 55. 

LEVEL II INPUT 

LEVEL 2 REOUIRES THE ITEMS REQUIRED FOR LEVEL 1 PLUS ALL OTHER 
ITEMS EXCEPT LOCATIONS 64, G5,' 91-100, AND 175-192. THE MINIMUM 
INPUT FQR THIS LEVEL WOULD BE LEVEL 1 DATA PLUS DAMPING, FLASH 
TEMPERATURE DATA, SOLUTION ITERATION DATA, AND PHASING CONSTANTS. 
DEFAULT VALUES ARE ZERO UNLESS INDICATED OTHERWISE REI.OW. 


LEVEL III INPUT 

LEVEL 3 REOUIRES ALL ITEMS TO BE INPUT, UNLESS DEFAULT VALUES ARE 
INDICATED. THE MINIMUM INPUT FOR THIS LEVEL WOULD BE LEVEL 2 
DATA . 


NOTE: 


FOR HIGH CONTACT RATIO GEARS (CR > 2.), THE USER MUST IN- 
PUT THE CR (I.OC 54 AND 55). PROGRAM CALCULATES RATIO 
IF CR < 2. 


** G E A R SYSTEM PARAMETERS ****************** 


LOC 

NAME 

* A 

* * * * 

1 

LEVEL 

2 

DP 

3 

PSANC 

4 

PSANGB 

5 

PS 10 

C 

N 1 

7 

N2 

S 

» o 
!«•> 

9 

FW 1 

10 

FW2 

11 

FW:J 

12 

N 

13 

i; 


DESCRIPTION 

********* 


******* 


TRIGGER FOR LEVEL OF INPUT DATA 
DIAMETRAL PITCH ( NORMAL PLANE) 

SUN-PLANET DRIVE SIDE PRESSURE ANCLE G PD 

RING-PLANET COAST SIDE PRESSURE ANGLE 6 PD 

SUN -PLANET COAST SIDE PRESSURE ANCLE G PD 

RING-PLANET DRIVE SIDE PRESSURE ANGLE £ PD 

** 0.0 IE PRESSURE ANGLES FOR DRIVE & COAST SIDES EQUAL 

HELIX ANGLE £ PD 

** 0.0 FOR SPUR GEARS 

NUMBER OF TEETH ON SUN GEAR 

PLANET GEARS 
RING GEAR 

AXIAL FACE WIDTH OF SUN GEAR 

PLANET GEARS 
RING GEAR 
OF PLANET GEARS 

PLANETARY SYSTEM, I .E. .RING GEAR FIXED 
STAR SYSTEM, I. E. .PLANET CARRIER FIXED 
NON PLANETARY, I. E. , NO RING AND NO CARRIER 

SUN- INPUT, PLANET-OUTPUT 
NON PLANETARY, I. E. , NO SUN AND NO CARRIER, 

PLANET INPUT, RING OUTPUT 
PLANETARY SYSTEM WITH FLEXIBLE PLANET CARRIER 
STAR SYSTEM WITH FLEXIBLE RING GEAR RIM 


NU* 

JBER 

IF 

K 1 

IF 

K 2 

* IF 

K ’ J 

IF 

K-5 

* IF 

KT, 

* IF 

K 7 
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******* K - C AND K - 7 CURRENTLY NOT AVAILABLE * * x * * * * A 


* IE K-G OR K-7 LOC S9 IS REQUIRED FOR PLANET CARRIER 
OR RING GEAR RIM STIFFNESS VALUE ALONG THE LINK— OI 7 — 
ACTION. 

** FOR K-4 AND RUNOUT ERROR. INPUT LOC 200. 


14 TORO 

15 RPM 

16 RPMF 

17 INTVL 


***] 

* GEAR 

22 

El 

23 

E2 

24 

E3 

25 

MU1 

2G 

MU2 

27 

M13 

* * * ! 

ft GEAR 

28 

MS 

20 

MC 

30 

MR 


AXIAL INPUT TORQUE ON SUN GEAR (OR PLANET IF K-D) 
INITIAL AXIAL ROTATIONAL SPEED OF SUN GEAR 

(OR PLANET IF K-5) 

FINAL AXIAL ROTATIONAL SPEED OF SUN GEAR FOR SPD RANGE 
NUMBER OF MAIN INTERVALS THE SPEED RANGE DIVIDED INTO 
** FOR ONE RPM ONLY , RPM- RPMF AND INTVL- 1.** 

STEP SIZE -(RPMF - RPM)/ INTVL. WILL BE DIVIDED 
BY 5 IN AREAS OF PEAK LOADS 

MATERIAL PROPERTIES ******************************* 

YOUNGS MODULUS * E-OC OF SUN GEAR (DEFAULT - 30.) 

PLANET GEARS “ « “ 

RING GEAR “ “ “ 

POISSONS RATIO OF SUN GEAR (DEFAULT - .30) 

PLANET GEARS “ “ “ 

RING GEAR “ “ “ 


EQUIVALENT MASSES ********************************* 


EQUIVALENT MASS OF SUN CEAR ABOUT ROTATIONAL AXIS 
PLANET CARRIER 
RING GEAR 


31-50 MP< I) 


PLANET GEARS 


*** NOTE: 


FOR K-G OR K-7.CAP.riEn OR RING MASSFOR TOTAL 
UNIT, NOT SEGMENTS. 


******* TRIGGER FOR DOUBLE HELICAL GEARING ********************* 

51 DB1IEL IF - 0.0 SINCLE HELICAL GEARING IF PSIO .NE. 0.0 
IF - 1.0 DOUBLE HELICAL GEARING IF PSIO .NE. 0.0 


* 


* 


******** 


****************** 


END OF LEVEL 
******* 


INPUT (UNLESS HIGH CONTACT RATIO) 
****************** 
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* " * 






******* GEAR meshing damping ratios **************************** 

52 ZSP DAMPING RATIO OF SUN-PLANET MESHES 

53 ZRP “ “ “ RING-PLANET MESHES 

* * * NOTE: IF CONTACT RATIO > 2, CRSP AND CRRP MUST DE INPUT, BUT 

IF ZERO IS INPUT THE PROGRAM WILL CALCULATE FOR CR < 2 ONLY. 

54 CRSP INVOLUTE CONTACT RATIO OF SUN-PLANET MESHES 

55 CRRP 44 “ 44 RING-PLANET MESHES 

******* ACTIVE FACE WIDTH VALUES ******************************* 
**NOTF: LOC 50 -. 59 FOR HELICAL GEARS ONLY 

5G LISP INACTIVE SUN-PLANET FACE WIDTH ON LEFT 

57 L2SP “ “ 44 “ “ “ RIGHT 

58 L1RP INACTIVE RING-PLANET FACE WIDTH ON LEFT 

59 L2RP “ “ “ “ “ “ RIGHT 

******* SUN-PLANET PROFILE MODIFICATION INPUT DATA ************* 

GO PCTSOD SD AS A PERCENT CF SOD (7.) 

61 PCTSOK SE AS A PERCENT OF SOE (7.) 

62 DELD DISENGAGEMENT TIP RELIEF <IN.), MINIMUM 

63 DELE ENGAGEMENT TIP RELIEF ( IN. >, MINIMUM 

**NOTE : SICN CONVENTION— POSITIVE INPUT, MATERIAL REMOVED 

****** XNVDX AND XNVEX ARE NOT REPUIRED FOR LEVEL 2 * * * * 

64 XNVDX DISENGAGEMENT PROFILE MODIFICATION SHAPE FACTOR 

65 XNVEX ENGAGEMENT 

********************************* 

GG PMODSP PERCENT OF TIP RELIEF ON Sl'N GEAR, ENGAGEMENT 

67 PMDSPD PERCENT OF TIP RELIEF ON SUN GEAR, DISENGAGEMENT 

- 7. TIP RELIEF ON PLANET 

68 DLSP TOLERANCE AT START OF PROFILE MODIFICATION 

78 DLTOL TOTAL TIP RELIEF TOLERANCE (SUN * PLANET 

AND/OR RING * PLANET, DEPENDING ON SYSTEM TYPE) 

SOD - LENGTH OF DISENGAGEMENT PART OF THE LINE OF ACTION 
SOE - 44 44 ENGAGEMENT 44 44 44 44 44 44 

SD & SE ARE THE SEGMENTS OF THE SOD £ SOE THAT ARE USED, WHEN 
PROFILE MODIFICATIONS ARE MADE. 

******* RING-PLA.NET PROFILE MODIFICATION INPUT DATA ************ 

69 PTS0D3 SD AS A PERCENT OF SOD 

70 PT.S0K3 SE AS A PERCENT OF SOE 

71 DELD3 DISENGAGEMENT TIP RELIEF 

72 DEI.E3 ENGAGEMENT TIP RELIEF 
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i-r- 40#* 
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* * * * * LOCATION 73 AND 74 NOT REQUIRED FOR LEVEL 2 * * * * * 


- 7 . TIP RELIEF ON RING GEAR 
TOLERANCE AT START OF PROFILE MODIFICATION • 
TOTAL TIP RELIEF TOLERANCE, FOR BOTH SUN ♦PLANET 
AND/OR RING ♦ PLANET 

****** MES1I MODIFICATION DUE TO FACE WIDTH CROWNING *********** 


“ DISENGAGEMENT FOR SP MES! 


73 

XNVDX3 

74 

XNVEX3 

* * * 

* * * * 

75 

PMODRP 

76 

PMDRPD 

77 

DLRP 

78 

DLTOL 

****** 

MESH Ml 

80 

LECSP 

81 

LDCSP 

82 

DLECSP 

83 ’ 

DLDCSP 

84 

LECRP 

85 

LDCRP 

86 

DLECRP 

87 

DLDCRP 

****** 

FLEX IB! 

88 

ZCR 

89 

CARI5K 

****** 

TOOTH ! 

on 

W I OK 

* * * 

* * 91 

01 

SPKSP 

92 

COMASP 

93 

COMBS P 

94 

COMCSP 

95 

COMDSP 

96 

SPKRP 

97 

COM ARP 


M “ DISENGAGEMENT FOR RP MESI 


IF K-C, PLANET CARRIER DAMPING RATIO 
IF K-7, RING CEAR RIM DAMPING RATIO 

IF K-G PLANET CARRIER STIFFNESS A LON G - THE - 1. 1 N E - OF - 
ACTION (///IN) 

IF K-7 RING CEAR RIM STIFFNESS ALONC-TIIE-LINE-OF- 
ACTION (/'/IN) 


IF - 0.0 PREPROCESSOR CALULATES TOOTH PAIR COMPLIANCE 
IF - 1.0 INPUT COMPLIANCE DATA IN LOCATIONS 01-101 

* * * * * 91 - 100 NOT REQUIRED FOR LEVEL 2************ 

SINGLE TOOTH PAIR SPR INGRATE OF Sl'N -PLANET MESHES 
COMPLIANCE CONSTANT <S/Sf!>**l OF SUN- PLANET MESHES 

2 

3 

4 

SINCLE TOOTH PAIR SPRINCRATE OF RING-PLANET MESHES 
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9R COMURP 
99 COMCRP 
100 COMDRP 


2 

3 

4 


WHERE S - IS THE LENGTH ALONC THE LINE OF ACTION TO THE POINT 
OF ENGAGEMENT 

AND SO - THE LENGTH OF THE LINE OF ACTION 


ft ft ft 

« « * * 

* « 

101 

HRTZSP 

HE! 

102 

HRTZP.P 



** THE 

PRO! 


PLIANC1 


STRESS 

10S 

XPSPSP 

IF 



IF 

104 

XPSPRP 

IF 



IF 

111 

CONVEK 

IF 



IF 

******* 

TOOTH PAIR 

120 

BYPASS 

IF 



IF 

121 

RADIOS 

MAX 

122 

RADTOP 


123 

RADTOK 


124 

RADTIS 

MAX 

125 

RAD7IP 


12G 

RADTIR 


127 

TOOTHS 

NOM 

128 

TOOTH!' 


129 

TOOTHR 


120 

RADFIS 

MAX 

131 

RADFIP 


122 

RADFIR 



RING-PLANET 


STRESS 

STRAIN 


AXIAL SEGMENTS 
■ 1.0 CONVECTIVE COMPLIA' 
IN LOCATIONS 701 TO SOO 


RING- PLANET 
RING-PLANET 


121-122 


PLANET GEARS 
RING CEAR 


PLANET GEARS 
RINC GEAR 


PLANET GEARS 
RING GEAR 


( ROOT 
RADII'S) 


PLANET GEARS 
RINC GEAR 
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****** 

FLASH 

TEMPERATURE DATA ********************************* 

140 

SROU 

SURFACE ROUGHNESS (RMS, MICRO INCHES) 

141 

TOIL 

OIL INLET TEMPERATURE 

142 

CFMAT 

MATERIAL CONSTANT 

143 

OILTYP 

TYPE OF OIL IN CEARBOX 


1 FOR MIL-L-23690 

2 FOR MIL-L- 7S0S 

3 TO BE DETERMINED 

4 TO BE DETERMINED 

(0 DEFAULTS TO TYPE 1) 

******* SOLUTION ITERATION DATA ******************************** 


100 TOLER ITERATION TOLERANCE FOR BOUNDARY CONDITION CONVERGENCE 

(7. * 100, DEFAULT FOR LEVEL 1 - .01) 

101 NLOOP NUMBER OF ITERATIONS FOR BOUNDARY CONDITION CONVERGENCE 

NLOOR MAXIMUM - 20, DEFAULT - 20. 

****** TOOTH PAIR GEOMETRIC TOLERANCE DATA ******************** 


ICO 

1C1 

1C2 

1G3 

164 

165 
ICC 
1G7 


DLROS TOLERANCE ON TIP RADIUS OF SUN CEAR (DEFAULT - .002) 

DLROP PLANET GEARS 

DLROR RING CEAR “ “ “ 

BSE EDGE BREAK ON TOPLAND OF ALL GEARS (DEFAULT - .010) 

DLRR ROOT RADIUS TOLERANCE OF ALL GEARS (DEFAULT - .005) 

DLCDN CENTER DISTANCE TOLERANCE (TOWARD TIGHT MESH) ALL GEARS 
DLCDP CENTER DISTANCE TOLERANCE (AWAY FROM TIGHT MESH) 

DLMBT MACHINE BACKLASH TOLERANCE ALI. GEARS (DEFAULT - .002) 


NOTE: THE TOOTH ROOT FILLET RADIUS TOLERANCE 


.005 IN THE PREPROCESSOR 


******* 3-DIME.NSIONAL FACTORS FOR SPUR GEARS ONLY ************** 


* * * * * 175 - 192 NOT REQUIRED FOR LEVEL 2************ 


175 

176 

177 

178 

179 

180 


SPRNRS 


SPRNBP 


SPRNBR 


SPRNLS 

SPRNI.P 

SPRNLR 


IF - 0. RIM BENDING EFFECTS NOT INCLUDED IN SUN CEAR 
(DEFAULT - 10. E* 10) 

IF > 0. BENDINC SPRISCRATEt IN-///RAD ) OF SUN RIM 
IF - 0. RIM BENDINC EFFECTS NOT INCLUDED IN PLANET CEAR 
(DEFAULT - 10 . E* 10 ) 

IF > 0. BENDING SPRINCRATE( IN-///RAD) OF PLANET RIM 
IF " (J. RIM BENDING EFFECTS NOT INCLUDED IN RING GEAR 
(DEFAULT’- 10.EM0) 

IF > 0. BENDING SPRINCRATFM IN-///RAD) OF RING RIM 

RADIUS FROM RIM TO PD (REQUIRED WITH SPENDS) SUN 

RADIUS FROM KIM TO I’D (REQUIRED WITH SPRNBP) PLANET 

RADIUS FROM RIM TO RD (REQUIRED WITH SPRNBR) PLANET 


181 

182 

183 

181 

I8r, 

18fi 

187 

183 

lS'J 

IPO 

191 

192 


EFWP1S 


EFV.D2S 

EFVD1P 


EFV.D2P 


FI- V D 1 1! 


FFV.D2R 

STRSFS 

STRSCS 

STRSEP 


STRSCP 

STKSER 

STRSCR 


EFFECT I VF SIN FACE WIDTH FACTOR AT TOOTH TIP 
IF - 0 . ( DEFAULT - 1.0) 

EFFECTIVE SIN FACE WIDTH FACTOR AT FII.LET 
IF - U. (DEFAULT - 1.0) 

EFFECTIVE PLANET FACE WIDTH FACTOR AT TOOTH TIP 
IF - 0. ( DEFAULT - 1.0) 

EFFECTIVE PLANET FACE WIDTH FACTOR AT FILLET 
IF - 0. (DEFAULT - 1.0) 

EFFECTIVE KING FACE WIDTH FACTOR AT TOOTH TIP 
IF - 0. (DEFAULT - 1.0) 

EFFECTIVE RING FACE WIDTH FACTOR AT FII.LET 
IF - 0. (DEFAULT - 1.0) 


STRESS PISTRIL . FACTOR FOR END EFFTS 
IF - 0. (DEFAULT - 1.0) 

STRESS DISTRID. FACTOR FOR END EFFTS 
IF - 0. (DEFAULT - 1.0) 

STRESS D I STRIP. FACTOR FOR END EFFTS 
IF - 0. (DEFAULT - 1.0) 


AT SUN WIDTH EDGE 


AT SIN WIDTH CENT 
AT PLAN WDTil EDGE 


STRESS DISTRIR. FACTOR FOR 
IF - 0. (DEFAULT - 1.0) 
STRESS DISTRIR. FACTOR FOR 
IF - 0. ( DEFAULT - 1.0) 
STRESS D I STRIP. FACTOR FOR 


END 

EFFTS 

AT 

PLAN 

WDTii 

CENT 

END 

EFFTS 

AT 

RING 

V.DTH 

EDGE 

END 

EFFTS 

AT 

RING 

U DTI! 

CENT 


IF - 0. (DEFAULT - 1.0) 


****** * 


HELIX ANGLE ERRORS ************************************* 


195 

DELl’SP 

HELIX ANGLE ERROR FOR SIN -PLANET MESH 


190 

DEI.PPP 

“ " “ “ RING-PLANET MESH 


197 

DELPS2 

IF Dl*H!:i. - 0.0 DELPS2 - 0.0 
IF DllHEL - l.U DEI.PS2 - HELIX ANGLE ERROR OF 
HALF OF SUN-PLANET DOIT.LE HELICAL 

RIGHT 

GEARS 

193 

DEI.PR2 

IF DDIIFI. • 0.0 DELPR2 - 0.0 

IF DPHFL - 1.0 DFI.PP2 - HELIX ANCLE ERROR OF 

RIGHT 


HALF OF RING-PLANET DOUBLE HELICAL GEARS 
******* TOOH! PAIR SPACING ERRORS ****************************** 


20(1 DR SI n RUNOUT ERROR FOR EXTERNAL-EXTERNAL SINGLE MESH ONLY 
(K*4), DISPLACEMENT ERROR OF SUN CENTER 

*>NOTI: FOR D0». r.LJ: HEI.ini. GEARS ( DP.HEL - 1.0) THE FIRST 10 

VALUES OF I IN THE E(I,J) ARRAYS ARE FOR THE LEFT HALT 
01 THE DOCILE HELICAL GEARS A‘D THE LAST 10 LOCATIONS 
ARC FOR THE RIGHT HALF OFTIIF DOUBLE HELICAL GEARS, 

WHERE I - THE PLANET NIMHER. 
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ORIGINAL PAGl 

OF POOR QUALITY 


221 TO 210 ESP( 1 , 1 )SCN-PLANET ERROR ARRAY FOR TOOTH 1 


2*1 1 TO 2C0 ESP<I,2) 2 
2G1 TO 280 ESP( 1,3) 3 
231 TO 300 ESP< 1,4) 4 
301 TO 320 ESP( 1,5) 5 

321 TO 340 ERP( 1,1) RING-PLANET ERROR ARRAY FOR TOOTH 1 
341 TO 3G0 EKP( 1,2) 2 
3C1 TO 380 ERP< 1,3) 3 
381 TO 400 ERP( 1,4) 4 
401 TO 420 ERP( 1,5) 6 


******* PLANET GEARS PHASING CONSTANTS' ************************* 

421 TO 410 KSP(I) SlIN-PLANET PHASING CONSTANTS ARRAY 

441 TO 4C0 KRP< I ) RING-PLANET PHASING CONSTANTS ARRAY 

******* INITIAL BOUNDARY CONDITIONS **************************** 


481 TO 500 
G01 TO 520 
521 TO 540 
541 TO 5G0 


XOSP(I) SIN-PLANET INITIAL DISPLACEMENT BOUNDARY COND 
Xur.P(I) RING-PLANET INITIAL DISPLACEMENT BOUNDARY COND 
XISP(I) Sl.N-PLANET INITIAL VELOCITY BOUNDARY COND 
XIRP(I) RING-PLANET INITIAL VELOCITY BOUNDARY COND 


******* PROGRAM PLOT SELECTIONS FOR SPUR GEARS ONLY ************ 


G01 

PLTLD 

IF 

- 

0 0 




IF 

- 

1.0 

LOAD PLOTS, NORMALIZED LOAD — DYNAMIC/STATIC 

602 

PLTPV 

IF 

- 

0.0 




IF 

- 

1.0 

PV PLOTS, PRESSURE SLIDING VELOCITY 

GOO 

PLTHS 

IF 

- 

0.0 




IF 

- 

1.0 

HFRT7. STRESS PLOTS 

G04 

PLTI'T 

IF 

- 

o.n 




IF 

- 

1.0 

FLASH TEMPERATURE PLOTS 

GOG 

PLTSS 

IP 

- 

0.0 




IF 

- 

1.0 

SIN GEAR ICYV.OOD STRESS PLOTS 

GOC 

PLTPS'S 

IF 

- 

0.0 




IF 

- 

1.0 

SUN -PLANET PLANET GEAR HE WOOD STRESS PLOTS 

G07 

PLTPRS 

IF 

- 

0.0 




IF 

- 

1.0 

R INC- PLANET PLANET GEAR HEW. 001) STRESS PLOTS 

COS 

PLTRPS 

IF 

- 

o.n 




IF 

- 

l.n 

RING CHAP HEYROOD STRESS PLOTS 


79 


«ct!rr 


****** * PROGRAM CHECK RUN i ************************************** 


099 Clirri; IF - 0.0 REGULAR Rt.N TO COMPLETION 

IF - 1.0 EXITS PROGRAM BEFORE DYNAMIC SOLUTION 

TO ALLOW CHECKING OF INPUT DATA AND THE 
PREPROCESSOR RESULTS 

******* CONVECTIVE COMPLIANCE MATRIX FOR HELICAL GEARING ******* 

** CURRENTLY NOT AVAILABLE WITHOUT PROGRAM ALTERATIONS, IE. COMMENT 
STATEMENTS MUST RE REMOVED 


7U1 

TO 

709 

LI.-1 

TO 

9 

L- 1 

711 

TO 

719 

LL-1 

TO 

9 

L-2 

721 

TO 

729 

LL-1 

TO 

9 

L-3 

721 

TO 

720 

LL-1 

TO 

9 

L-4 

711 

TO 

749 

LL-1 

TO 

9 

L-5 

751 

TO 

759 

LL-1 

TO 

9 

L-G 

7GI 

TO 

7G9 

LL-1 

TO 

9 

L-7 

771 

TO 

779 

LL-1 

TO 

9 

L-8 

7SI 

TO 

729 

LL-1 

TO 

9 

L-9 

* * 

* 1 

» * * 

* * 

* * 

* * 

* * * 


J 


J 


***************** 


***** OUTPUT DESCRIPTION 

THE FOLLOWING LIST INDICATES THE INFORMATION AND RESULTS THAT 
MAY APPEAR IN THE OITPLT, IN THE ORDER THEY WILL APPEAR. 


CASE TITLE— THE TITLE AND/OR OTHER INFORMATION IM'tT ON Till! FIRST 
DATA CARD. 


TIP MODIFICATION— IF THERE IS INSUFFICIENT TIP CLEARANCE. AN 
ADJUSTMENT IS MADE INTERNALLY AND A MESS AGE OITI'IT UNTIL 
SUFFICIENT CLEARANCE II OBTAINED. 

INVOLUTE MODIFICATION TABLE — SHOWS TUI. PROCESSED RESULTS OF 
ANY INPUT PROFILE MODIFICATIONS. TABLES ARE PRINTED FOR 
ENGAGEMENT AND DISENGAGEMENT WHICH INCLUDE MINIMIM AND 
MAXIMUM INVOLUTE MODIFICATIONS , DIAMETER AND CORRESPONDING 
ROLL ANCLE FOR THE MODIFIED PORTION OF THE PROFILE. THIS 
WILL BE OUTPUT FOR THE SI'S -PLANET MESH AND/OR RING-PLANET 
MESH IN CONJUNCTION WITH THE CORRESPONDING INPUT DATA. 

INPUT DATA— THE INPUT DATA AND PRLPROCENSED GEOMETRIC DATA IS 
PRINTED FOR SIN-PLANET MILS II AND/OR RING-PLANE! MESil. 


00 





amCNNAL PAGE tS 
OF POOR QUALITY 


FLEXIBILITY— IF THE FLEXIBLE PLANET CARRIER OPTION IS IN EFFECT 
A MESSAGE APPEARS THAT INDICATES THIS. 

ADDITIONAL INPIT DATA— NUMBER OF PLANETS, EQUIVALENT MASSES, ETC. 

COMPLIANCE CONSTANTS— CALCULATED CONSTANTS FOR FOURTH ORDER 
COMPLIANCE EQUATION. 

ITERATED BOUNDARY CONDITIONS— BOUNDARY CONDITIONS PRINTED, FOL- 
LOWED BY THE CURRENT SPEED BEING EXAMINED AND CORRESPONDING 
MAXIMUM LOADS (FOR SUN-PLANET AND/OR RING-PLANET MESHES). 

IF A SPEED SURVEY WAS RUN, THE SPEED CORRESPONDING TO THE 
OVERALL MAXIMUM LOAD FROM THE RANGE CALCULATED IS OUTPUT 
WITH THE MAXIMUM LOAD. THIS IS THE SPEED USED FOR THE 
REMAINING CALCULATIONS, I.E. STRESS. 

MAXIMUM VALUES— TABLE(S) OF MAXIMUM VALUES CALCULATED ARE PRINTED 
FOR EACH SUN-PLANET MESH AND/OR RING-PLANET MESH. THESE 
WILL APPEAR FOR THE NO ERROR SOLUTION AND EACH ERROR SOLU- 
TION. TWO ABBREVIATIONS. APPEAR— PV-PRESSURE SLIDING VELO- 
CITY AND PD-PITCH DIAMETER. 

ERROR MATRIX— IF TOOTH SPACING ERRORS ARE INPUT, OR GENERATED 
FOR RUNOUT SOLUTION, A TABLE OF ERRORS IS PRINTED BEFORE 
THE TABLES OF MAXIMUM VALUES. 
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APPENDIX B: EXAMPLE PROBLEM AND DATA SET 


EXAMPLE PROBLEM 


Diametral Pitch * 8.4667 
Pressure Angle * 22.5° 

Number of Sun Teeth * 14 

Number of Ring Teeth » 70 

Number of Planet Teeth « 28 

Face Width, Sun & Planets = 1.1811 in. 

Face Width, Ring * 1.4184 in. 

Input Torque * 30 in. lb. 

Input RPM = 7000 rpm 



Planet Gear 
Ring Gear, fixed 
Sun Gear 
Planet Carrier 
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Example Data Set 
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CALCULATES DYNAMIC LOAD RESPONSE FOR 

A. SINGLE MESH SPUR GEARING 

B. MULTIPLE MESH SPUR GEARING. 

STAR OR PLANETARY 

C. SINGLE MESH HELICAL GEARING 




APPENDIX C: EXAMPLE OUTPUT 
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THIS PROGRAM HAS DEVELOPED AT HAMILTON STANDARD DIVISION OF UNITED TECHNOLOGIES , WINDSOR LOCKS 
J. PIKE. R. CORNELL. W. WESTERVELT, AND L. EOlO. FUNDING HAS GIVEN UNDER CONTRACT BY NASA-LE’.JIS 
CLEVELAND AND N.DNI TOPED DY D. TCNMSEKO. 


NOT ENOUGH TIP CLEARANCE 0.0047 CHANGE RA1P BY .005 
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ROTATIONAL PLANE 
INVOLUTE MODIFICATIONS 
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y. 


NO. TEETH - S'JN 14.0000 NUTBER OF TEETH 14.0000 28.0000 


UKKjWAL FAGiL IS 

OF POOR QUALITY 
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NUMBER OF PLANETS 3 

NUMBER OF BOUNDARY CONDITION ITERATIONS 20 

TOLEPANCE FOR BCUIIDART CONDITION CONVERGENCE 0.10000E-01 



EQUIVALENT MASS OF SUN GEAR 0.280S3E-01 

EQUIVALENT HSC3 GF PLANE I CARRIER 0.l'i459E»02 

EQUIVALENT NASS OF RING GEAR 0.0 

EQUIVALENT NASS OF PLANET » 1 0.77535E-01 
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RTHG- PLANET ; MESHES 1 THRU N, LEFT TO RIGHT 

HEPATIC l» E?.'0 DISPLACEMENT 13 0.572<i3E-0* 0.S1126E-05 #.<»090SE-05 

SV’TUlS DISPLACEMENT 0.6<M34E-06 0.5131iE-05 0.9C552E-0S 

ITCPATIC-I. E!3l:.“J VELOCITY 13 -0.2:t3t«-02 -0.16«S7E-01 0.K.171E-01 

st;pt:?:3 velccitt -o.c^i76E-o: -o.i?o3?e-oi o.i*c9;c-oi 


mVIMUM lO'P FC 3 SL’H-P4.*::£t = 19.13 
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maximum values for sun-planet mesh 
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THE EFFECTIVE COMTACT PATIO = 1.4500 
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